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Abstract 

Thermodynamical properties of nuclear matter undergoing multifragmentation are 
studied within a simplified version of the statistical model. An exact analytical solution 
has been found for the grand canonical ensemble. Excluded volume effects are taken into 
account in the thermodynamically self-consistent way. In thermodynamic limit the model 
exhibits a first order liquid-gas phase transition with specific mixed phase properties. 
An extension of the model including the Fisher's term is also studied. The possibility 
of the second order phase transition at or above the critical point is demonstrated. The 
fragment mass distributions in the different regions of the phase diagram are discussed. 
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I. Introduction 



Nuclear multifragmentation, i.e. multiple production of intermediate mass fragments, is 
one of the most spectacular phenomena in intermediate energy nuclear reactions. The statis- 
tical multifragmentation model (SMM) has been developed during last two decades (see ffl, |2| 
and references therein). Numerous comparisons with experimental data (see e.g. Ref. @. [|) 
show that it is rather successful in explaining many important features of nuclear multifrag- 
mentation. Moreover, there are serious indications j|, [5| that multifragmentation in equi- 
librated systems is related to a liquid-gas phase transition in nuclear matter. Recently a 
simplified version of the SMM has been proposed || [7| to study this relationship. The calcu- 
lations within the canonical ensemble of non-interacting fragments suggest the existence of a 
1-st order phase transition. This was demonstrated by increasing the total number of particles 
in the system up to A=2800. These numerical calculations appeared to be rather efficient due 
to the recursive formula || which explicitly expressed the canonical partition function for A 
particles in terms of the canonical partition functions for 1,2,...,A-1 particles. The recursive 
procedure essentially simplifies the numerical evaluations and makes them possible for rather 
large values of A. However, the analytical studies of the system behavior in the thermody- 
namic limit, i.e. when both the system mass number A and volume V go to infinity, are still 
missing. On the other hand, the investigation of the thermodynamic limit is crucial to proof 
the existence of a phase transition and to study its exact nature. 

In the present paper we give an exact analytical solution of the simplified version of the 
SMM performing its complete study in the thermodynamic limit. The accurate treatment of 
the excluded volume effects is an important part of our study. We work in the grand canonical 
ensemble which significantly simplifies the mathematical problems and leads to the explicit 
analytical solution in the thermodynamic limit. 

The paper is organized as follows. In Sect. II the simplified version of the SMM is 
formulated. In Sect. Ill we introduce the isobaric partition function. The relationship between 
its singularities and the phase transition existence is discussed in Sect. IV. Sect. V is dealing 
with the first order phase transition and properties of the mixed phase. In Sect. VI the 
possibility of a second order phase transition is demonstrated. Fragment mass distributions 
are discussed in Sect. VII. Sect. VIII is reserved for conclusions. 



II. Model Formulation 

Let us consider a system composed of many different nuclear fragments and characterized 
by the total number of nucleons A, volume V and temperature T. The system states are 
specified by the multiplicities {n k } {n k = 0,1,2,...) of /c-nucleon fragments {k = 1,2,..., A). 
The partition function of a single fragment with k nucleons is assumed to have the standard 
form (here and below the units h = c = 1 are used) 

fmTk\ 3/2 

UJk = V l~2 — J Zk ' ^ 

where m is the nucleon mass and the fragment mass is approximated asmfc. The first two 
factors in Eq. ([[]) originate from the non-relativistic thermal motion of the fc-fragment in the 
volume V at temperature T. The last factor, z k , represents the intrinsic partition function of 
the /c-fragment. For k = 1 (nucleon) we take Z\ — 4 (4 internal spin-isospin states) and for 
fragments with k > 1 we use the expression motivated by the liquid drop model (see details 



in Ref. §): 

z k = exp ( - ^ 1 , (2) 




where 

fk = -[Wo + T 2 /e D ] k + a(T)0 3 + r Tln(fe) (3) 

is the internal free energy of the fc-fragment. Here W D = 16 MeV is the volume binding energy 
per nucleon, T 2 /e Q is the contribution of the excited states taken in the Fermi-gas approxima- 
tion (e Q = 16 MeV) and cr(T) is the surface free energy tension which is parameterized in the 
following form 

/T 2 — T 2 \ 5 ^ 

°F) = -o [ff-^) m-n (4) 

with cr = 18 MeV and T c = 18 MeV. Finally, tT ]n(k) is the phenomenological Fisher's term 
(r is a dimensionless constant) which we introduce to generalize our discussion. This form 
of z k (with r = 0) was used in Refs. |5[ 0. The symmetry and Coulomb contributions to the 
free energy are neglected. Such a model, however, appears to be a good starting point for the 
phase transition studies. 

III. Isobaric Partition Function 

The canonical partition function (CPF) of the ideal gas of nuclear fragments has the 
following form: 



Za(V,T) = E II ZTJ S(A-J2kn k ) , (5) 

{n k } k k ' k 

where ^-function takes care of the total baryon number conservation. An important assump- 
tion of the model is that the fragments do not overlap in a coordinate space. This gives rise 
to the repulsive interaction which we take into account in the Van der Waals excluded volume 
approximation. This is achieved by substituting the total volume V in Eq. @ by the free 
volume Vf = V — bA, where b = l/p Q (p Q = 0.16 fm~ 3 is the normal nuclear density). The cal- 
culation of the CPF d5|) is difficult because of the restriction J2k ^ n k — A. This constraint can 
be avoided by calculating the grand canonical partition function (GCPF). Using the standard 
definition we can write 

Z{V,T,n) = £exp(^) Z A (V-bA,T)e(V - bA) (6) 
= E II ^7 \( V - bJ2 kn k) MT) expQjik/T)] 9 ( V - bJ2kn k ) , 

{n k } k nk - L\ k ) J \ k ) 

where <j> h (T) = u k (T, V)/V. 

The presence of the theta-function in the GCPF (^J) guarantees that only configurations 
with positive value of the free volume are counted. However, similarly to the delta function 
restriction in Eq. (H), it makes again the calculation of Z(V,T,fi) ([]) to be rather difficult. 
This problem can be solved |1(J by performing the Laplace transformation of Z(V, T, p). This 
introduces the so-called isobaric partition function (IPF): 
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After changing the integration variable V —>■ V, the constraint of 0-function has disappeared. 
Then all were summed independently leading to the exponential function. Now the inte- 
gration over V in Eq. (0) can be straightforwardly done resulting in 



where 
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Here we have introduced the shifted chemical potential v = // + W + T 2 /e Q . Now the 
GCPT can be obtained by the inverse Laplace transformation 



I ra+ioo 

Z(y,T,fi) = — ds exp(sV) Z(s,T,fi) 

LH% J a—ioc 



(10) 



where the integration is taken along the imaginary axis. In the thermodynamic limit V — > oo 
the system pressure is defined by the farthest-right singularity, s*(T, /i), of IPF Z(s,T,fi) (0) 
(see Ref. [IT] for details): 



V— >oo 1/ 



The simple connection of the farthest-right s-singularity of Z, Eq. (8), to the asymptotic, 
V — > oo, behavior of Z, Eq. (6), is a general mathematical property of the Laplace transform. 
Due to this property the study of the system behavior in the thermodynamic limit V — > oo 
can be reduced to the investigation of the singularities of Z. 

IV. Singularities of IPF and Phase Transitions 

The IPF, Eq. (0), has two types of singularities: 
1). the simple pole singularity defined by the equation 



s g (T,fj,) = F{s g ,T,n) , 



(12) 



2). the essential singularity of the function J-(s, T, fi) it-self at the point si where the coefficient 
in linear over k terms in the exponent is equal to zero, 



v 

Tb 



(13) 



The simple pole singularity corresponds to the gaseous phase where pressure is determined 
from the following transcendental equation 



. /mT\ 3/ 2 

p ' (T,Al) = UrJ T 



Z\ exp 
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(14) 



The singularity si(T,p) of the function jF(s, T, p) @ defines the liquid pressure 



7^ 

6 



(15) 



In the considered model the liquid phase is represented by an infinite fragment, i.e. it 
corresponds to the macroscopic population of the single mode k = oo. Here one can see the 
analogy with the Bose condensation where the macroscopic population of a single mode occurs 
in the momentum space. 

In the (T, /i) -regions where v < bp g (T,p) the gas phase dominates (p g > pi), while the 
liquid phase corresponds to v > bp g (T, p). The liquid-gas phase transition occurs when two 
singularities coincide, i.e. s g (T,p) = sj(T, p). A schematic view of singular points is shown in 
Fig. la for T < T c , i.e. when a > 0. The two-phase coexistence region is therefore defined by 
the equation 

Pi(T,p) = p g (T,fi) , i.e., v = bp g (T,p). (16) 

One can easily see that ^(s, T, fi) is monotonously decreasing function of s. The necessary 
condition for the phase transition is that this function remains finite in its singular point 
si = u/Tb: 

F{s h T,n) < oo . (17) 

The convergence of F is determined by r and a. At r = the condition flTTD requires cr(T) > 0. 
Otherwise, !F(si, T, p) = oo and the simple pole singularity s g (T, p) ( |i~2D is always the farthest- 
right s-singularity of Z (0) (see Fig. lb). At T > T c , where u(T) = 0, the considered system 
can exist only in the one-phase state. It will be shown below that for r > 5/2 the condition 
(|T7j) can be satisfied even at cr(T) = 0. 



V. 1-st Order Phase Transition and Mixed Phase 

At T < T c the system undergoes the 1-st order phase transition across the line p* 
defined by Eq. (|i~6|) . Its explicit form is given by the expression (W = W Q + T 2 /e Q ): 



p*(T) 
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The points on the line p*(T) correspond to the mixed phase states. We first consider the case 
when r = 0. The line p*(T) ([18]) for this case is shown in Fig. 2. 

The baryonic density is calculated as (dp/dp)T and is given by the following formulae in 
the liquid and gas phases, respectively: 
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where the function p^ is defined as 
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(21) 



Due to the condition flTB]) this expression is simplified in the mixed phase: 
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(22) 

This formula clearly shows that the bulk (free) energy acts in favor of the composite fragments, 
but the surface term favors single nucleons. 

Since at a > the sum in Eq. ( |2"2"D converges at any r, p id is finite and according to Eq. ( PDD 
p g < 1/b. Therefore, the baryonic density has a discontinuity Ap = pi — p g > across the 
line p*{T) ([18|) shown for r = in upper panel of Fig. 2. The discontinuities take place also 
for the energy and entropy densities. The phase diagram of the system in the (T, p)-plane is 
shown in lower panel of Fig. 2. The line p*{T) (0) corresponding to the mixed phase states is 
transformed into the finite region in the (T, p)-plane. In this mixed phase region of the phase 
diagram the baryonic density p is a superposition of the liquid and gas baryonic densities: 



P 



A pi + (1 - A) Pg . 



(23) 



Here A (0 < A < 1) is a fraction of the system volume occupied by the liquid inside the mixed 
phase. Similar linear combination is also valid for the energy density: 



with {i — l,g): 
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The pressure on the phase transition line p*(T) (|18|) is a monotonously increasing function 
of T 



p*(T)= P g(T,p*(T)) 
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(28) 



Fig. 3 shows the pressure isotherms as functions of the reduced density p/p D for r = 0. Inside 
the mixed phase the obtained pressure isotherms are horizontal straight lines in accordance 
with the Gibbs criterion. These straight lines go up to infinity when T — > T c — 0. This formally 
corresponds to the critical point, T = T c , p = p c = 1/b and p c = oo, in the considered case of 



r = 0. For T > T c the pressure isotherms never enter into the mixed phase region. Note that, 
if cx(T) would never vanish, the mixed phase would extend up to infinite temperatures. 

Inside the mixed phase at constant density p the parameter A is temperature dependent 
as shown in Fig. 4: A(T) drops to zero in the narrow vicinity of the boundary separating 
the mixed phase and the pure gaseous phase. This specific behavior of A(T) causes a strong 
increase of the energy density (|24"D and as its consequence a narrow peak of the specific heat 
Cy(T) = (de/dT)p. It should be emphasized that the energy density is continuous at the 
boundary of the mixed phase and the gas phase, and the sharpness of the peak in Cy is 
entirely due to the strong temperature dependence of A(T) near this boundary. A narrow 
peak of the specific heat was observed in the canonical ensemble calculations of Refs. || 0. 
However, in contrast to the expectation in Refs. || [7|, the height of the Cy{T) peak is not 
equal to infinity and its width is not zero in the thermodynamic limit considered in our study. 
Note also that the shape of the Cy(T) depends strongly on the parameter r and on the chosen 
value of the baryon density p. 

VI. Possibility of 2-nd Order Phase Transition 

The results presented in Figs. 2-4 are obtained for r = 0. New possibilities appear at 
non-zero values of the parameter r. At0<r<5/2 the qualitative picture remains the same 
as discussed above, although there are some quantitative changes. For r > 5/2 the condition 
(|T7|) is also satisfied at T > T c where <r(T) = 0. Therefore, the liquid-gas phase transition 
extends now to all temperatures. Its properties are, however, different for r > 7/2 and for 
r < 7/2 (see Fig. 5). If r > 7/2 the gas density is always lower than 1/6 as p^ is finite. 
Therefore, the liquid-gas transition at T > T c remains the 1-st order phase transition with 
discontinuities of baryonic density, entropy and energy densities. The pressure isotherms as 
functions of the reduced density p/p Q are shown for this case in Fig. 6. 

At5/2<r<7/2 the baryonic density of the gas in the mixed phase, Eqs. ( pC||2"2"D , becomes 
equal to that of the liquid at T > T c , since p id — > oo and p™ %x = 1/b = p Q . It is easy to prove 
that the entropy and energy densities for the liquid and gas phases are also equal to each other. 
There are discontinuities only in the derivatives of these densities over T and p, i.e., p(T,p) 
has discontinuities of its second derivatives. Therefore, the liquid-gas transition at T > T c for 
5/2 < t < 7/2 becomes the 2-nd order phase transition. According to standard definition, the 
point T = T c , p = 1/b separating the first and second order transitions is the tricritical point. 
One can see that this point is now at a finite pressure. Fig. 7 shows the pressure isotherms as 
functions of the reduced density p/p Q . 

It is interesting to note that at r > the mixed phase boundary shown in Fig. 5 is not so 
steep function of T as in the case r = presented in Fig.3. Therefore, the peak in the specific 
heat discussed above becomes less pronounced. 

VII. Fragment Mass Distributions 

The density of fragments with k nucleons can be obtained by differentiating the gas pressure 
(|14}) with respect to the fc-fragment chemical potential pk — kp. This leads to the fragment 
mass distribution P(k) in the gas phase 

P g {k) = a k 3/2 ~ T exp(- ai k - a 2 k 2/3 ) , (29) 

where ai = (bp g — v)jT > 0, a 2 = a /T and a is the normalization constant. Since the 
coefficients a ,ai,ci2 depend on T and p the distribution P(k) (|29|) has different shapes in 



different points of the phase diagram. In the mixed phase the condition ( pp leads to a\ — 
and Eq. (^) is transformed into 



P™™(k) = a Q k 3/2 ~ T exp(-a 2 fc 2/3 ) . (30) 



The liquid inside the mixed phase is one infinite fragment which occupies a fraction A of 
the total system volume. Therefore, in a large system with A nucleons in volume V (A/V = p) 
the mixed phase consists of one big fragment with XV p Q nucleons (liquid) and (1 — X)Vp g nu- 
cleons distributed in different fc-fragments according to Eq. (|30|) (gas) . At low T most nucleons 
are inside one big liquid-fragment with only few small gas-fragments distributed according to 
Eq. (|30D with large a 2 . At increasing temperature the fraction of the gas fragments increases 
and their mass distribution becomes broader since a 2 (T) in Eq. ( PD| ) decreases. Outside the 
mixed phase region the liquid disappears and the fragment mass distribution acquires an ex- 
ponential falloff, Eq.(p9l). Therefore, the fragment mass distribution is widest at the boundary 
of the mixed phase. At even higher temperatures, T > T c , the coefficient a 2 vanishes. 

Details of the fragment mass distribution depend on the parameter r. At r < 5/2 we 
observe a sudden transformation of the large liquid fragment into light and intermediate mass 
fragments in the narrow vicinity of the mixed phase boundary. This sudden change of the 
fragment composition has the same origin as a narrow peak in the specific heat, i.e. a sharp 
drop of A(T) near the mixed phase boundary (see Fig. 4). For larger r all these changes are 
getting smoother. 

An interesting possibility opens when 5/2 < r < 7/2. As shown in Fig. 7 the mixed phase 
in this case ends at the tricritical point T = T c , p = p . In this point both the coefficients ai 
and a 2 vanish and the mass distribution becomes a pure power law 



3 



P g (k) = a Q k*- r . (31) 
At t > 7/2 the mixed phase exists at all T. Thus the mass distribution of gaseous fragments 



inside the mixed phase fulfills a power-law (|31|)at all T > T c . 
VIII. Conclusions 

We have used a simplified version of the statistical multifragmentation model (SMM) |T| 
to establish the relationship between multifragmentation phenomenon and a liquid-gas phase 
transition in nuclear matter. Recently, in Refs. H [7) interesting peculiarities of this model were 
found numerically in the canonical ensemble formulation. In present paper this simplified SMM 
is solved analytically by considering the thermodynamic limit V — > 00 in the grand canonical 
ensemble. The progress has been achieved by applying an elegant mathematical method which 
reduces the description of phase transitions to the investigation of singularities of the isobaric 
partition function. In this way we have exactly solved the model in the thermodynamic limit. 
The excluded volume effects are fully taken into account. 

The model clearly demonstrates the 1-st order phase transition of the liquid-gas type. It 
is rather surprising that in thermodynamic limit the liquid phase emerges as an infinite-mass 
fragment. The structure of the mixed phase and some peculiar properties near its boundary 
are discussed in details. The phase diagram appears to be rather sensitive to the value of 
the parameter r in the Fisher's free energy term included in our treatment. New interesting 
possibilities for the phase diagram emerge for r > 5/2 in comparison with the case when 
r < 5/2. The case 5/2<r<7/2is particularly interesting because of the appearance of the 
tricritical point separating the 1-st and 2-nd order phase transitions. 



The results presented in this paper will be further developed taking into account additional 
physical inputs (e.g. finite size effects, Coulomb interactions and symmetry energy) to make 
the model closer to reality. 
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Fig. 1. Schematic view of singular points of the Isobaric Partition Function, Eq. (8), at 
T < T c (a) and T > T c (b). Full lines show F{s,T, fx) as a function of s at fixed T and 
/i, /ii < /i2 < A*3 < A*4- Dots and asterisks indicate the simple poles (s g ) and the essential 
singularity of function T it-self (s/). At /X3 = fJ>*(T) the two singular points coincide signaling 
a phase transition. 




Fig. 2. Phase diagram in T — p (upper panel) and T — p (lower panel) planes for r = 0. 
The mixed phase is represented by the line p*(T) in upper panel and by the extended region 
in the lower panel. Liquid phase (shown by crosses) exists at density p = p Q . Point C is the 
critical point. 




Fig. 3. Pressure isotherms (thin solid lines) as functions of reduced density p/p for r = 0. 
The isotherms are shown for T = 4, 8, 10, 14, 18 and 22 MeV from bottom to top. The bound- 
ary of the mixed and gaseous phases is shown by the thick solid line. Liquid phase is indicated 
by crosses. The critical point T = T c = 18 MeV, p = p Q corresponds to infinite pressure. 
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Fig. 4. Volume fraction of the liquid phase A(T) for r = shown as a function of temperature 
at densities p/p Q = 1/6, 1/3, 1/2, 2/3, 5/6 (from bottom to top). 
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Fig. 5. Phase diagrams in T — p plane for r = 3.6 (upper panel) and r = 2.6 (lower panel). 
Point C in the lower panel is the tricritical point. Crosses correspond to the liquid phase of 
the first order phase transition and dots correspond to the states of the second order one. 




Fig. 7. The same as Fig. 3, but for r = 2.6. Point C in the upper panel is the tricritical 
point. The lower panel differs by the scale. Crosses correspond to the liquid phase of the first 
order phase transition and dots correspond to the states of the second order one. 



